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ABSTRACT 

Near-infrared (NIR) spectra of the subluminous Type la supernova SN 1999by are 
presented which cover the time evolution from about 4 days before to 2 weeks after 
maximum light. Analysis of these data was accomplished through the construction of 
an extended set of delayed detonation (DD) models covering the entire range of normal 
to subluminous SNe la. The explosion, light curves (LC), and the time evolution of 
the synthetic spectra were calculated self-consistently for each model with the only free 
parameters being the initial structure of the white dwarf (WD) and the description of 
the nuclear burning front during the explosion. From these, one model was selected 
for SN 1999by by matching the synthetic and observed optical light curves, principly 
the rapid brightness decline. DD models require a minimum amount of burning during 
the deflagration phase which implies a lower limit for the ^^Ni mass of about O.IMq 
and consequently a lower limit for the SN brightness. The models which best match 
the optical light curve of SN 1999by were those with a ^^Ni production close to this 
theoretical minimum. The data are consistent with little or no interstellar reddening 
[E{B—V) < 0.12”^) and we derive a distance or 11±2.5 Mpc for SN 1999by in agreement 
with other estimates. 

Without any modification, the synthetic spectra from this subluminous model match 
reasonably well the observed IR spectra taken on May 6, May 10, May 16 and May 24, 
1999. These dates correspond roughly to —4 d, 0 d, and 6 d and 14 d after maximum 
light. Prior to maximum, the NIR spectra of SN 1999by are dominated by products 
of explosive carbon burning (O, Mg), and Si. Spectra taken after maximum light are 
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dominated by products of incomplete Si burning. Unlike the behavior of normal Type la 
SNe, lines from iron-group elements only begin to show up in our last spectrum taken 
about two weeks after maximum light. The implied distribution of elements in velocity 
space agrees well with the DD model predictions for a subluminous SN la. Regardless 
of the explosion model, the long duration of the phases dominated by layers of explosive 
carbon and oxygen burning argues that SN 1999by was the explosion of a white dwarf 
at or near the Chandrasekhar mass. 

The good agreement between the observations and the models without fine-tuning 
a large number of free parameters suggests that DD models are a good description of 
at least subluminous Type la SNe. Pure deflagration scenarios or mergers are unlikely 
and helium-triggered explosions can be ruled out. However, problems for DD models 
still remain, as the data seem to be at odds with recent 3-D models of the deflagration 
phase which predict significant mixing of the inner layers of the white dwarf prior to 
detonation. Possible solutions include the effects of rapid rotation on the propagation 
of nuclear flames during the explosive phase of burning, or extensive burning of carbon 
just prior to the runaway. 

Subject headings: supernovae - individual: SN1999by 


1. Introduction 

On April 30, SN 1999by was independently discovered at about IS™' by Arbor et al. (1999), 
and Papenkova, Li, & Filippenko (1999). Images with the Katzman Automatic Imaging Telescope 
(Treffers et al. 1997; Li et al. 1999) of the same field provide an upper limit of 19.3™ on April 25.2 
UT. The supernova (SN) was found in the Sb galaxy NGC 2841 which has been the host of three 
previous supernovae; SN 1912A, SN 1957A, and SN 1972R (Papenkova, Li, & Filippenko 1999). 
Based on optical spectra, SN 1999by was identified as a Type la SN (Gerardy & Fesen 1999). 
Shortly thereafter, Garnavich, Jha, & Kirshner (1999) reported that SN 1999by showed stronger 
than normal Si II 5800 A absorption and depressed flux near 4000 A, suggesting that SN 1999by 
would be a significantly subluminous Type la event. 

According to Bonanos et al. (1999), SN 1999by reached a maximum light of niB = 13.80™ ± 
0.02™ on UT 1999 May 10.5, and a maximum in the V band of my = 13.36™ ± 0.02™. The 
Lyon/Meudon Extragalactic Database (LEDA)^, gives the heliocentric radial velocity of NGC 2841, 
corrected for Virgo infall, as 811.545 km s“^. Using ffg = 65 km s“^ Mpc“^ puts the distance to 
NGC 2841 at 12.5 Mpc with a distance modulus of 30.5™. Using this distance modulus and the 
photometry of Bonanos et al. (1999), the absolute peak magnitudes of SN 1999by are Mb = 
— 16.68™ and My = —17.12™. Thus, SN 1999by is underluminous by roughly 2.5 magnitudes as 
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compared to a typical Type la supernova (SN la). 

Detailed analysis of the optical light curves (LCs) by Toth & Szabo (2000) confirm these 
basic results. The light curve of SN 1999by shows a very steep post-maximum brightness decline 
My {AM /^t=i5d) of 1.35”^ to 1.45™. Based on detailed fits of the LCs, Toth & Szabo (2000) find 
the interstellar reddening to be E{B — V)< 0.1™ which is in agreement with the values for galactic 
reddening given by Schlegel et al. (1998, E{B — V) = 0.015™), and Burstein & Heiles (1982, 
E{B — F) « 0™). Recently, Garnavich et al. (2001) provided detailed optical LCs and spectra and 
found values consistent with previous data. 

These measurements imply that SN 1999by was as underluminous as SN 1991bg, the proto¬ 
typical example of the subluminous Type la subclass (Filippenko et al. 1992; Leibundgut et al. 
1993). Other members include SN 1992K (Hamuy et al. 1996a,b), SN 1997cn (Turatto et al. 1998), 
and SN 1998de (Modjaz et al. 2000). Some defining characteristics of the subclass are rapidly 
declining light curves (MB(AM^t=i 5 rf) ~ 1.9™), peak magnitudes fainter than normal by 2-3 mag, 
and redder colors at maximum light {{B — F)max — 0.4-0.5™). One of the currently active areas 
in SNe la research concerns the nature of these subluminous events. Theoretical interpretations 
of subluminous SN la include all three types of explosion mechanism: the centrally triggered det¬ 
onation of a sub-Chandrasekhar mass WD, deflagration and delayed detonations of massive WD 
and two merging WDs (see sect. 3 &: Woosley &: Weaver, 1994 ; Hoflich et ah, 1995; Nugent et ah, 
1995, Milne et ah 1999). The possibility has also been raised that SN1991bg-like SNe la should 
not be classified with other SNe la at all as they may arise from different progenitors. 

SN 1999by is one of the best observed SNe la with data superior to that of SN 1991bg. In 
addition to the studies of optical light curves and spectra mentioned above, detailed polarization 
spectra of SN 1999by have been obtained and analyzed (Howell et ah 2001). Whereas ‘normal’ 
SNe la tend to show little or no polarization (e.g. Wang, Wheeler, &: Hoflich 1997), this supernova 
was significantly polarized, up to 0.7%, indicating an overall asphericity of w 20%. This result 
suggests that there may be a connection between the observed asphericity and the subluminosity 
in SNe la. 

In recent years it has become apparent that infrared spectroscopic observations can be used as 
valuable tools to determine the chemical structure of SN la envelopes. For instance, near-infrared 
spectra can be used to locate the boundaries between explosive carbon and oxygen burning, or 
between complete and incomplete silicon burning by measuring Mg II, Ca II and iron-group lines 
(Wheeler et ah 1998). However, the difficulties involved in obtaining high quality infrared spectra of 
supernovae has limited IR studies to a small handful of objects. The study of SN 1994D by Meikle 
et ah (1996) is the only one to include both pre- and post-maximum spectra. Since most of their 
spectra was of the 1-1.3 //m region, Meikle et ah (1996) concentrated on determining the origin of 
a feature at 1.05 /xm. They concluded that this feature might be due to either He I 1.083 /xm or 
Mg II 1.0926 ^m, but found difficulties with both identifications. Detailed modeling of SN 1994D 
and SN 1986G was able to reproduce the basic NIR spectral features and their evolution (Hoflich 
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1997; Wheeler et al. 1998). In particular, the identification of the 1.05 //m feature in SN 1994D as 
due to Mg II could be established (recently confirmed by Lentz et al. 2001), and a broad feature, 
appearing between ~ 1.5 and 1.9 /um was identified as a blend of iron-group elements. 

We note that the nature of subluminous SNe la is important for the use of SNe la as distance 
indicators. In particular, the calibration of the brightness decline relation depends critically on 
subluminous SNe la because they significantly increase the required baseline observations. From 
this perspective, it may turn out to be critical to determine whether all SNe la form a homogeneous 
class of objects or not. However, current limitations of model calculations will not allow us to 
improve the accuracy of the current estimates for the absolute distance to SN 1999by, and the 
possible implications of our results on the use of SNe la for cosmology is not the subject of this 
paper. 

In this paper, we present near-infrared spectra of SN 1999by covering the time evolution of 
the supernova from about four days before to two weeks after maximum light. Detailed models of 
the SN explosion based on a delayed detonation (DD) scenario are used to analyze the data. The 
explosion models, light curves, and synthetic spectra are calculated in a self-consistent manner. 
Given the initial structure of the progenitor and description of the nuclear burning front, the 
light curves and spectra are calculated from the explosion model without any further tuning. The 
observations and data reduction are discussed in §2. In §3, general properties of the explosion 
models, light curves and spectra are discussed for normal and subluminous SNe la within the DD 
scenario. In §4, a specific model for this SN is chosen by matching the properties of the predicted 
light curve to those observed in optical studies of SN 1999by. Synthetic optical and IR spectra for 
this model are then presented and the latter are compared to the observed IR spectra. Detailed 
line identifications for the infrared features are provided. The effects on the predicted spectra of 
large scale mixing during the deflagration phase also examined. In §5, we discuss the implications 
for progenitors and alternative explosion mechanisms. We close in §6 with a final discussion of the 
results for SN 1999by, and put our findings in context with other SNe la and different explosion 
scenarios. 


2. Observations 

Low-dispersion (R ~ 700), near-infrared spectra of SN 1999by from 1.0-2.4 /xm were obtained 
with the 2.4m Hiltner Telescope at MDM Observatory during the nights of 6-10 May 1999. The data 
were collected using TIFKAM (a.k.a. ONIS), a high-throughput infrared imager and spectrograph 
with an ALLADIN 512 x 1024 InSb detector. This instrument can be operated with standard J, H, 
and K filters for broadband imaging, or with a variety of grisms, blocking filters, and an east-west 
oriented Ol'G slit for low {R ^ 700) and moderate {R ^ 1400) resolution spectroscopy. 

SN 1999by was observed with three different spectroscopic setups, which covered the 0.96- 
1.80 /xm, 1.2-2.2 ^m, and 1.95-2.4 ^m wavelength regions. The observations were broken into 
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multiple 600 s exposures. Between each exposure, the supernova light was dithered along the slit 
to minimize the effect of detector defects and provide first-order background subtraction. Total 
exposure times vary from spectrum to spectrum and are listed in the log of NIR observations in 
Table 1. 

Wavelength calibration of the spectra was achieved by observing neon, argon, and xenon arc 
lamps. The spectra were corrected for telluric absorption by observing A stars and early G dwarfs 
at similar airmass, chosen from the Bright Star Catalog (Hoffleit &: Jaschek 1982). Applying the 
procedure described by Hanson, Reike, & Luhman (1998), stellar features were removed from the 
G dwarf spectra by dividing by a normalized solar spectrum (Livingston & Wallace 1991; Wallace, 
Hinkle, & Livingston 1993)^. The results were used to correct for telluric absorption in the A stars. 
Hydrogen features in the corrected A star spectra were removed from the raw A star spectra and 
the results were used to correct the target data for telluric absorption. [For further discussion of 
this procedure see Hanson, Reike, &: Luhman 1998; Hanson, Conti, & Rieke 1996, and references 
therein.] The instrumental response was calibrated by matching the continuum of the A star telluric 
standards to the stellar atmosphere models of Kurucz (1993). 

The first three nights (6-8 May) were photometric, and J, H, and K broadband images of 
SN 1999by and Persson et al. (1998) photometric standards. The resulting photometry is listed in 
Table 1 and was used to set the flux levels of the corresponding spectra. A square bandpass was 
assumed and the flux levels thus attained are believed accurate to ~ 20-30%. For May 9 and 10 no 
flux information was available, but the observed spectral bands had large overlaps, and the relative 
flux levels were set by matching the data in the overlapping regions. 

Near-infrared spectra of a star 291^5 north and lT.% west of SN 1999by (Papenkova, Li, & 
Filippenko 1999) were also obtained and reduced in the same manner. This star was then used as a 
local relative spectrophotometric standard to reduce spectroscopic data from later observing runs. 
Spectra of SN 1999by covering 1.0-1.8 /xm were collected on five more nights during the following 
two weeks (16 May, 18 May, 20 May, 21 May and 24 May). The first three were obtained by P. 
Martini and A. Steed using TIFKAM on the 2.4m telescope at MDM, and the last two by S. Sakai 
using TIFKAM on the 2.1m telescope at the Kitt Peak National Observatory. 

Figure 1 shows a plot of all the NIR spectra. Since absolute fluxes were only available for 
the first three spectra, the data are presented in arbitrary flux units, and they have been shifted 
vertically for clarity. Regions of very low S/N due to strong telluric absorption have been omitted. 
The epochs listed are relative to 10 May 1999, the date of VMax given by Bonanos et al. (1999). 
Like all the subsequent plots, the data have been corrected to the 638 km s“^ redshift of NGC 2841 
(de Vaucouleurs et al. 1991). 


®NSO/Kitt Peak FTS data used here were produced by NSF/NOAO. 
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3. Models for the Explosion, Light Curves, and Spectra 

There is general agreement that SNe la result from some process involving the combustion 
of a degenerate white dwarf (Hoyle & Fowler 1960). Within this general picture, three classes of 
models have been considered: (1) An explosion of a CO-WD, with mass close to the Chandrasekhar 
limit, which accretes mass through Roche-lobe overflow from an evolved companion star (Whelan 
& Iben 1973). The explosion is triggered by compressional heating near the WD center. (2) An 
explosion of a rotating configuration formed from the merging of two low-mass WDs, caused by the 
loss of angular momentum through gravitational radiation (Webbink 1984; Iben & Tutukov 1984; 
Paczyhski 1985). (3) Explosion of a low mass CO-WD triggered by the detonation of a helium 
layer (Nomoto 1980; Woosley, Weaver, & Taam 1980; Woosley & Weaver 1986). 

Only the first two models appear to be viable for normal Type la SNe as the third, the sub- 
Chandrasekhar WD model, has been ruled out on the basis of predicted light curves and spectra 
(Hoflich et al. 1996b; Nugent et al. 1997). Still, theoretical interpretations of subluminous SN la 
include all three types of explosion mechanism (Woosley & Weaver 1994; Hoflich, Khokhlov, & 
Wheeler 1995; Nugent et al. 1995; Milne, The, & Leising 1999). The possibility has also been 
raised that the subluminous SN 1991bg-like supernovae may arise from different progenitors and 
should not be classified with other SNe la at all. 

Within Mch models, it is believed that the burning front starts as a subsonic deflagration. 
However, the time evolution of the burning front is still an open question. That is, whether 
the deflagration front burns through the entire WD (Nomoto, Thielemann, & Yokoi 1984), or 
alternatively, transitions into a supersonic detonation mode as suggested in the delayed detonation 
(DD) model (Khokhlov 1991; Woosley & Weaver 1994; Yamaoka et al. 1992). DD models have been 
found to reproduce the optical and infrared light curves and spectra of ‘typical’ SNe la reasonably 
well (Hoflich 1995, H95 hereafter; Hoflich & Khokhlov 1996; Hoflich et al. 1996b; Fisher et al. 1998; 
Nugent et al. 1997; Wheeler et al. 1998; Lentz et al. 2001). 

In addition, DD models provide a natural explanation for the brightness decline relation 
M{AM^t=i 5 d) observed in the light curves of SNe la (Phillips et al. 1987; Hamuy et al. 1995, 
1996b; Suntzeff et al. 1999). This is a consequence of the temperature dependence of the opacity 
and the lack of a dependence of the explosion energy on the amount of produced (Hoflich, 
Muller, & Khokhlov 1993; Khokhlov, Muller, &: Hoflich 1993; Hoflich et al. 1996b; Mazzali et al. 
2001). Thus, within the delayed detonation scenario, both normal and subluminous SNe la can be 
explained as variants of a single phenomenon (H95, Hoflich, Khokhlov, & Wheeler 1995, Umeda 
et al. 1999). This makes the DD model an attractive scenario and it is used as the basis for our 
present analysis. 

In our models, the explosions LCs and spectra are all calculated self-consistently. Given the 
initial structure of the progenitor and a description of the nuclear burning front, the light curves 
and spectra are calculated directly from the explosion model without any additional parameters. 
Since the predicted observables follow directly from the supernova model, this approach provides 
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a direct link between the observations and the explosion physics and progenitor properties. 

Finally, we note that detailed analyses of observed SN la spectra and light curves indicate that 
mergers (Benz et al. 1990) and pure deflagration SNe (such as the W7 model) could contribute to 
the population of bright SNe la (Hoflich & Khokhlov 1996, hereafter HK96; Hatano et al. 2000). 
On the other hand, it is not obvious that subluminous SNe la and the brightness decline relation 
can be understood within pure deflagration models. In particular, the amount of burning, and 
consequently the explosion energy, will be strongly correlated with production. Since no 

detailed studies for deflagration models have been performed, these models cannot be ruled out 
but, as we discuss in §5, the IR spectra of SN 1999by are hard to reconcile with a pure deflagration. 


3.1. Numerical Methods 

3.1.1. The Explosion 

Despite recent progress in our understanding of nuclear burning fronts, current 3-D models 
are not sufficiently evolved to allow for a consistent treatment of the burning front throughout all 
phases (Khokhlov 2001, see also §4.3 & §5). Thus detailed models rely on parameterized descriptions 
guided by 3-D results. Within spherical models, the transition from deflagration to detonation can 
be conveniently parameterized by a density ptr which has been found to be the dominating factor 
for the determining chemical structure, light curves and spectra. 

We have calculated explosion models using a one-dimensional radiation-hydro code (HK96) 
that solves the hydro dynamical equations explicitly by the piecewise parabolic method (Collela 
& Woodward 1984). Nuclear burning is taken into account using an extended network of 606 
isotopes from n,p to (Thielemann, Nomoto, & Hashimoto 1996, and references therein). The 

propagation of the nuclear burning front is given by the velocity of sound behind the burning front in 
the case of a detonation wave, and in a parameterized form during the deflagration phase calibrated 
by detailed 3-D calculations (e.g. Khokhlov 2001). We use the parameterization as described in 
Dominguez &: Hoflich (2000). For a deflagration front at distance rh^m from the center, we assume 
that the burning velocity is given by Ubum = 'fnax{yt,vi)., where vi and vt are the laminar and 
turbulent velocities with 

Vt = 0.2 y/ar g Lf, [ 1 ] 

with 

ar = (a — l)/(a + 1) 

and 

Ol = p (rburn)/P (^burn)- 

Here aT is the Atwood number, Tj is the characteristic length scale, and p'^ and p~ are the densities 
in front of and behind the burning front, respectively. The quantities a and Lj are directly taken 



from the hydrodynamical model at the location of the burning front and we take Lf = 'rbuj.n(t)- 
The transition density is treated as a free parameter. 


3.1.2. The Light Curves 

From these explosion models, the subsequent expansion, bolometric and broad band light 
curves (LC) are calculated following the method described by Hoflich, Wheeler, & Thielemann 
(1998), and references therein. The LC-code is the same used for the explosion except that 7 -ray 
transport is included via a Monte Carlo scheme and nuclear burning is neglected. In order to allow 
a more consistent treatment of the expansion, we solve the time-dependent, frequency-averaged 
radiation moment equations. The frequency-averaged variable Eddington factors and mean opac¬ 
ities are calculated by solving the frequency-dependent transport equations. About one thousand 
frequencies (in one hundred frequency groups) and about nine hundred depth points are used. At 
each time step, we use T(r) to determine the Eddington factors and mean opacities by solving the 
frequency-dependent radiation transport equation in a co-moving frame and integrate to obtain 
the frequency-averaged quantities. The averaged opacities have been calculated assuming local 
thermodynamic equilibrium (LTE). Both the monochromatic and mean opacities are calculated in 
the narrow line limit. Scattering, photon redistribution, and thermalization terms used in the light 
curve opacity calculations are taken into account. 

In previous works (e.g. Hoflich & Khokhlov (1996)), the photon redistribution and thermaliza- 
tion terms have been calibrated for a sample of spectra using the formalism of the equivalent two 
level approach (H95). Eor increased consistency, we use the same equations and atomic models but 
solve the rate equations simultaneously with the light curves calculation at about every 100 *^ time 
step, at the expense of some simplifications in the NLTE-part compared to H95. For the opacities 
we use the narrow line limit, and for the radiation fields we use the solution of the monochromatic 
radiation transport using w 1000 frequency groups. 


3.1.3. Speetral Calculations 

Our non-LTE code (H95, and references therein) solves the relativistic radiation transport 
equations in a co-moving frame. The spectra are computed for various epochs using the chemical, 
density, and luminosity structure and 7 -ray deposition resulting from the light curve coder. This 
provides a tight coupling between the explosion model and the radiative transfer. The effects of 
instantaneous energy deposition by 7 -rays, the stored energy (in the thermal bath and in ionization) 
and the energy loss due to the adiabatic expansion are taken into account. Bound-bound, bound- 
free and free-free opacities are included in the radiation transport which has been discretized by 
about 2 X 10^ frequencies and 97 radial points. 

The radiation transport equations are solved consistently with the statistical equations and 
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ionization due to 7 -radiation for the most important elements and ions. Typically, between 27 
to 137 bound levels are taken into account. We use C I (27/123/242), O 1 (43/129/431), Mg II 
(20/60/153), Si II (35/212/506), Ca II (41/195/742), Ti II (62/75/592), Fe II (137/3120/7293), 
Co II (84/1355/5396), Ni II (71/865/3064) where the first, second and third numbers in brackets 
denote the number of levels, bound-bound transitions, and number of discrete lines for the radiation 
transport. The latter number is larger because nearby levels within multiplets have been merged 
for the rates. 

The neighboring ionization stages have been approximated by simplified atomic models re¬ 
stricted to a few NLTE levels -I- LTE levels. The energy levels and cross sections of bound-bound 
transitions are taken from Kurucz (1993) starting at the ground state. The bound-free cross sections 
are taken from TOPBASE 0.5 as implemented in Strafiburg (Mendoza & Cunto 1993). Collisional 
transitions are treated in the ‘classical’ hydrogen-like approximation (Mihalas 1978) that relates 
the radiative to the collisional gf-values. All form factors are set to 1. 

About 10® additional lines are included (out of a line list of 4 x 10^) assuming LTE-level 
populations. The scattering, photon redistribution, and thermalization terms are computed with 
an equivalent-two-level formalism that is calibrated using NLTE models. 


3.2. Modeling Results 

Spherical dynamical explosions, light curves, and spectra are calculated for both normal bright 
and subluminous SNe la. Due to the one-dimensional nature of the models, the moment of the 
transition to a detonation is a free parameter. The moment of deflagration-to-detonation transition 
(DDT) is conveniently parameterized by introducing the transition density, ptr; at which DDT 
happens. 

Within the DD scenario, the free model parameters are: 1) The chemical structure of the 
exploding WD, 2) Its central density pc at the time of the explosion, 3) The description of the 
deflagration front, and 4) The density ptr at which the transition from deflagration to detonation 
occurs. Note that the density structure of a WD only weakly depends on the temperature because 
it is highly degenerate. The central density pc depends mainly on the history of accretion. Model 
parameters have been chosen for the progenitor and pc which allow us to reproduce light curves 
and spectra of ‘typical’ Type la supernovae. 

In all models, the structure of the exploding C/0 white dwarf is based on a model star with 5 
Mq at the main sequence and solar metallicity. Through accretion, this core has been grown close 
to the Chandrasekhar limit (see Model 5p0y23z22 in Dominguez, Hoflich, & Straniero 2001). At 
the time of the explosion of the WD, its central density is 2.0x10® g cm“^ and its mass is close 
to 1 . 37 M 0 . The transition density ptr has been identified as the main factor which determines the 
production and, thus the brightness of a SNe la (H95; Hoflich, Khokhlov, & Wheeler 1995; 
Iwamoto et al. 1999). The transition density ptr from deflagration to detonation is varied between 
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8 and 27 x 10® g cm“^ to span the entire range of brightness in SNe la. The models are identified 
by 5p0z22.ext where ext is the value of ptr in 10® g cm“^. Some of the basic quantities are given 
in Table 2. 


3.2.1. Explosion Models 

Here we will restrict our discussion to the basic features of DD models which are relevant for 
the spectral analysis of SN 1999by. For a more detailed discussion of the hydrodynamical evolution 
of DD-models, see Khokhlov (1991), H95, and references therein. 

During the deflagration phase, the WD is lifted, starts to expand, and after burning about 0.22 
to 0.28 Mq, the transition to a detonation is triggered. Figures 2 and 3 show the density, velocity 
and chemical structures for representative models after a homologous expansion has been achieved. 
The final velocity and density structures are rather similar. The expansion velocities decrease 
slightly with ptr because a significant amount of oxygen remains unburned and did not contribute 
to the energy production. In 5p0z22.8, oxygen remained unburned in the outer ~ O.4M0. Overall, 
the density is decreasing with radius because burning took place throughout the entire WD except 
in the very outer layers. In contrast, pure deflagration model such as W7 (Nomoto, Thielemann, 
& Yokoi 1984), pulsating delayed detonation (Hoflich, Khokhlov, & Wheeler 1995), and merger 
models (Khokhlov, Muller, & Hoflich 1993) show a significant density bump at the boundary 
between burned and unburned layers. 

Iron group elements are produced in the inner layers where density and temperature stay high 
for a sufficient time during the explosion, whereas intermediate mass elements are produced in the 
layers above where nuclear statistical equilibrium (NSE) has no time to set in. In the outer zones, 
products of explosive oxygen (e.g. Si, S) and explosive carbon burning (O, Mg, Ne) and some Si 
are seen. Only a very thin layer of unburned C/0 remains. 

For layers with final velocities v < 3000 km s“^, the densities in the early stages of the explosion 
are sufficiently high for electron capture. This results in the production of ®^’®®Fe & ®®Ni whereas 
in the outer layers only ®®Ni is produced. We employ the same description for the deflagration 
front in all models and, consequently, the very inner structures are almost identical in all models. 
A small dip in the Ni abundance is present at the DDT because the time scales for burning to NSE 
are comparable to the expansion time scales at densities of ~ 10^ g cm“^. This feature can be 
expected to be smeared out in multi-dimensional calculations (Livne 1999). 

After the DDT, a detonation front develops and the burning of fuel is triggered by compression 
of unburned material. Consequently, burning takes place under higher densities compared to a 
deflagration. In general, additional material is burned up to NSE and this extends the layers of 
®®Ni for all ptr > 1.2 x 10^ g cm“^. Eor models with smaller ptr, burning only proceeds up to Si 
despite the compression. This increases the production of intermediate mass elements (e.g. Si, S) 
at the expense of ®®Ni. In our models, the efficiency for ®®Ni production drops rapidly as a function 
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of decreasing ptr (see Table 2). We note that for ptr < Q x 10® g cm“^, nuclear burning stops all 
together for p < pth and, unlike any observed SN, one would see a very slowly expanding C/0 
envelope heated by little ®®Ni. 

There exists a rather strict lower limit for the absolute brightness of SNe la within the Chan¬ 
drasekhar mass models. This can be understood as follows. The absolute brightness of a SNe la 
is mainly determined by the ^^Ni. In general, is produced both during the deflagration and 
the detonation phase. For the most subluminous SNe la however, no is produced during the 
detonation. Within the Mch scenario and DD-models in particular, we need to burn a minimum 
of ~ O. 2 M 0 during the deflagration to achieve the required ‘lift/pre-expansion’ of the WD. During 
the deflagration phase, iron group elements are produced. The fraction of depends on the 
amount of material which undergoes electron capture, i.e. on the central density and, to a smaller 
extent, on the description of the burning front (Brachwitz 2001; Dominguez, Hoflich, & Straniero 
2001). For a model with pc = 2 x 10® g cm“^, we find a lower limit of Mi^^Ni) k, O.OOM©. Re¬ 
cently, Brachwitz (2001) studied the influence of pc between 1.5 to 8 x 10® g cm“® on the nuclear 
burning for similar DD-models with a high transition density which may describe ‘normal’ SN la. 
In those models, the produced during the deflagration phase varies by < 25% which provides 
an estimate for the possible variation in the minimum production. We conclude that there 
exists a lower limit of ^ O.IMq ± 25% within DD-models. 

In principle, 3-D effects may change the lower limit for the production. However, all 

current 3-D calculations show that rising blobs are formed at high densities which undergo complete 
Si-burning because the burning time scales to Ni shorter than the hydro time scales. In the case 
of a DDT or a transition to a phase of very fast deflagration, the amount of burning will likely 
be similar to the 1-D case. Thus, we do not expect a significant change in the minimum hQNi 
production. 


3.2.2. General properties of theoretieal light eurves 

Optical light curves and spectra of SN la are not the subject of this paper but we used their 
general properties for selecting a model for SN1999by. 

Within the DD scenario, both normal and subluminous SNe la can be produced (H95, Hoflich, 
Khokhlov, & Wheeler 1995, HKW95 hereafter; Umeda et al. 1999). The amount of ®®Ni (Msejvi) 
depends primarily on ptr (H95; HKW95; Umeda et al. 1999) and to a much lesser extent on the 
assumed value of the deflagration speed, initial central density of the WD, and initial chemical 
composition (ratio of carbon to oxygen). Models with smaller transition density give less nickel 
and hence both lower peak luminosity and lower temperatures (HKW95, Umeda et al. 1999). In 
DDs, almost the entire WD is burned and the total production of nuclear energy is almost constant. 
This and the dominance of ptr for determining the ^^Ni production are the basis of why, to first 
approximation, SNe la appear to be a one-parameter family. 
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The observed Mv'(AMAf=i 5 d) relation can be understood as an opacity effect (Hoflich et al. 
1996b), namely, the consequence of rapidly dropping opacities at low temperatures (Hoflich, Muller, 
& Khokhlov 1993; Khokhlov, Muller, &: Hoflich 1993; Mazzali et al. 2001). Less Ni means lower 
temperature and consequently, reduced mean opacities because the emissivity is shifted from the 
UV toward longer wavelengths with less line blocking. A more rapidly decreasing photosphere 
causes a faster release of the stored energy and, as a consequence, steeper declining LCs with 
decreasing brightness. The DD models give a natural and physically well-motivated origin of the 
My (AMAt=i 5 (i) relation of SNe la (Hamuy et al. 1995, 1996b; Phillips et al. 1987; Suntzeff et al. 
1999). 

The broad band B and V light curves are shown in Figure 4. Some general properties are 
given in Table 2 and Figure 5. We note that the decline of the B and V LC during the early 
tail is slightly smaller for the most subluminous models (e.g. 5p0z22.08) compared to the 
intermediate models (e.g. 5p0z22.16) because, at day 40, the y-rays are almost fully trapped in 
the former whereas, for the latter, the escape probability of y-rays is already increasing (HKW95). 
With varying ptr, we find rise times between 13.9 to 19.0 days, and the maximum brightness My 
spans a range —17.21™' and —19.35™ with B — V between -1-0.47™ to —0.03™. The value of My 
is primarily determined by the transition density. Variations in the progenitors (main sequence 
mass and metallicity) or the central density at the point of ignition causes an spread around My of 
Ri 0.2 ... 0.3™ (Hoflich, Wheeler, & Thielemann 1998; Dominguez, Hoflich, & Straniero 2001). As 
discussed above, the brightness decline relation My(AMAi=i 5 d) can be understood as an opacity 
effect or, more precisely, it is due to the drop of the opacity at low temperatures. Therefore, 
subluminous models are redder and the decline is steeper. In our series of models, My{/S.M^t=ibd) 
ranges between 1.45™ and 0.91™. However, My {AM/^t=i 5 d) is not a strictly linear relation (Fig. 
5). For bright SNe la (—19.35™ < My < —18.72™), called Branch-normals ((Branch, Romanishin, 
& Baron 1996; Branch 1999)), the linear decline is relatively flat, followed by a steeper decline 
toward the subluminous models. For normal-luminosity models, the photospheric temperature is 
well above the critical value at which the opacity drops (~ 7000 ...8000iL; Hoflich, Muller, & 
Khokhlov 1993) whereas the opacity drops rapidly in subluminous models soon after maximum 
light. 


4. Analysis of SN1999by 
4.1. Selection of the SN1999by Model 

Detailed B and V light curves of SN 1999by have been collected and published by Toth & Szabo 
(2000) and optical spectra around maximum light are available from observations at McDonald 
Observatory (Howell et al. 2001). Recently, SN 1999by spectra and LCs have been provided by 
Garnavich et al. (2001). Within our series of models, we selected the most suitable model based on 
the optical photometry. Unfortunately, there is a gap in the observed data starting about one week 
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after maximum light, a phase which is most suitable for selecting the appropriate models. Thus, 
as the main discriminator, we use the complete LCs as reconstructed by Toth & Szabo (2000). 

SN 1999by showed a decline rate My(AMAt=i 5 d) between 1.35™ and 1.45™ (Toth & Szabo 
2000) clearly ruling out all but the most subluminous models with transition densities of 8 and 
10 X 10® g cm“^ and maximum brightness of —17.21™ and —17.35™, respectively (Table 2). Both 
models are consistent with the early light curves given in Toth & Szabo (2000) i and Garnavich 
et al. (2001) (Fig.8). From the models and the early LC, the time of the explosion can be placed 
around April 27 ± 1 d, 1999. The model post-maximum LC in B is less reliable because it de¬ 
pends very sensitively on the size and time evolution of line blanketing and the thermalization. 

is 1.64™ and 1.73™ for 5p0z22.10 and 5p0z22.8, respectively. This compares rea¬ 
sonably well with Mb{^M/ s,t=i^d) = 1.87™ observed for SN1999by (Howell et al. 2001). In the 
LCs (Fig. 8), discrepancies of about 0.3™ show up after the time gap in the observations. At these 
times, the envelope of a subluminous SN la becomes transparent, and strong emission features 
emerge (see below). Thus these discrepancies are likely a consequence of discretization errors due 
to the use use of only 1000 frequency for the LC calculations. 

According to Bonanos et al. (1999), B — V = 0.44 ± 0.04™ near maximum light on May 10, 
1999. Garnavich et al. (2001) give B — V = 0.50 ± 0.03™. We hnd B — V = 0.47™ and 0.42™ for 
5p0z22.8 and 5p0z22.10, respectively. In previous studies, we hnd intrinsic uncertainties in B-V to 
be ~ 0.05™ for our models at maximum light (Hohich 1995). Thus, the interstellar reddening can 
be expected to be Eb-v < 0.11™ consistent with the values for our galaxy given by Burstein &: 
Heiles (1982) and Schlegel et al. (1998) and the estimates for the total extinction derived by Toth 
& Szabo (2000). For our models, if we assume my = 13.10 ± 0.05™, we then derive a distance 
modulus of My — my = 30.39™ ± 0.12™ (-0 to 0.35™). The hrst and second error-terms originate 
from the observational uncertainties plus the brightness range of the models and the reddening 
correction, respectively. ^From our models, we can derive a distance to NGC 2841 of 11 ± 2.5 Mpc 
and 12±1 Mpc if we include and neglect interstellar reddening, respectively. Both from observations 
(Schlegel et al. 1998) and our models, very small reddening is preferred. Following Schlegel et al. 
(1998), we adopt E{B — V) = 0.015™ for all comparisons with the observations. 

Based on the density, temperature and chemical structure predicted from the explosion model, 
'y — ray transport and light curve calculations, we have constructed detailed NLTE-spectra for 11 d, 
15 d, 22 d and 29 d after the explosion (days —4, 0, +7 and -|-14 after maximum light). For the 
subluminous model with the transition density ptr = 8 x 10® g cm“^ (Fig. 2), the evolution of the 
temperature, energy deposition by y-rays, and Rosseland optical depth is given as a function of 
radius (Fig. 6). Obviously, the total optical depth drops rapidly with time from 98, 40, 15 to 6 
due to the cooling; much faster than quadratic as would be expected from the geometric dilution 
alone (98, 52, 25 and 12), and thus more rapidly than in normal-bright SNe la (Hohich, Khokhlov, 
& Wheeler 1993; Khokhlov, Muller, & Hohich 1993)). As mentioned above, the reason is the 
rapid drop in the mean opacity for T < 6000...8000Ar which causes the outer layers to become 
almost transparent. Consequently, the photosphere recedes quickly in mass and even starts to 
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shrink in radius just a few days after the explosion. The corresponding expansion velocities at the 
photosphere are approximately 14000, 10500, 6500 and 4000 km s“^ at the four epochs computed. 

The low photospheric temperatures are due to the small amount of and its low expansion 
velocity which also causes the local trapping of 7 -rays (Fig. 6 ). In contrast to normal SNe la, most 
of the decay energy is deposited well below the photosphere up to about 2 weeks after maximum 
light. This explains, in part, the low ionization stages seen in the IR-spectra (see below). 

The location (and thus the velocity) of the photosphere determines which layers in the su¬ 
pernova envelope will form the spectral features. Before maximum light, the spectra sample the 
layers of explosive carbon burning which are O, Mg, Ne and Si-rich (Fig. 3). Up to 2 weeks after 
maximum light, spectra are formed within layers of incomplete silicon burning which are Si, S and 
Ca-rich. Only thereafter are the layers which are dominated by iron group elements finally exposed. 
This is very different from models for normal luminosity SNe la for which the iron rich layers are 
already exposed near maximum light. These profound differences are key for our understanding of 
the optical and IR spectra of subluminous SNe la and the differences compared to normal SNe la. 

Optical spectra are not the main subject of this paper, and detailed spectra have not been 
available for this analysis except for maximum light (Fig. 9, Howell et al. 2001). However, for 
completeness, the time sequence of optical spectra is given in Figure 7. Overall, the spectra show 
an evolution typical for SN la, (e.g. SN 1994D, H95). They can be understood as a consequence 
of the declining temperature and total optical depth with time, and the rapidly dropping temper¬ 
ature. The Doppler shift of the absorption component decreases rapidly because of the receding 
photosphere. For subluminous models, the entire envelope becomes almost transparent and the 
spectra start to enter the nebular phase at about two weeks after explosion. This marks the end of 
the applicability of our atomic models. Some features (see Fig. 9) are typical of subluminous SNe la 
and are a consequence of the low luminosity and photospheric temperature around maximum light. 
The Si II line at 5800 A is strong relative to the line at 6150 A, and the flux near 4000 A was 
depressed as in other subluminous SNe, and features due to Ti II were strong as was the O I line 
at 7500 A (Garnavich, Jha, & Kirshner 1999). As already suggested by Nugent et al. (1997), the 
strong Ti II absorption and the relative strength of the Si lines can be understood as a temperature 
effect. Within DD models with Mch progenitors, the strong O I line is a direct consequence of the 
fact that the spectra at maximum light are formed in massive layers of explosive carbon burning. 

While the structures and synthetic spectra of the models 5p0z22.8 and 5p0z22.10 are very 
similar, 5p0z22.8 is slightly preferred from a consideration of the expansion velocities at the pho¬ 
tosphere which are 10500 and 11500 km s“^, respectively. Therefore, we will will adopt 5p0z22.8 
in the following. Based on parameterized LTE model atmospheres (SYNOW, Fisher et al. 1999), a 
detailed analysis of the Doppler shifts of lines has been provided by Garnavich et al. (2001). They 
found Si II in the entire range between 11,300 down to 6500 km s“^ which is consistent with our 
hydrodynamical models (see Fig. 6 ). 
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4.2. Analysis of the IR-spectra of SN 1999by 

Near IR-spectra (0.9 to 2.5 ^m) of 5p0z22.8 at day 11, 15, 22 and 29 after explosion are 
compared to those of SN1999by at May 6, 10, 16, and 24 (Figs, refir.may6-13). The spectra are 
characterized by overlapping P-Cygni and absorption lines, and emission features. Unless noted 
otherwise, the observed wavelengths quoted the emission component. The overall energy distribu¬ 
tion in the continuum is determined by bound-free and Thomson opacities and the temperature 
evolution. In general, the spectra and their time evolution agree reasonably well with the model 
predictions and all the major features can be identified. 

Few detailed analyses and line identifications of IR spectra have been performed (Meikle et 
al. 1996; Bowers et al. 1997; Hoflich 1997; Wheeler et ah 1998; Hernandez et ah 2000) and these 
were limited to spectra of SN 1994D, SN 1998bu and post-maximum spectra of SN1986G. To put 
our analysis given below in context with normal-luminosity SNe la, we will refer to those works as 
reference points. 


4-2.1. Early time IR-spectra 

The observed NIR spectra from May 6 and 10 (Figs. 10 & 11) are quite similar, and comparison 
with the model spectra show that both of these spectra are formed in layers of explosive carbon 
burning. Spectral features can be attributed to C I, O I, Mg II and Si II. Typically, several thousand 
overlapping lines contribute to the overall opacity but few strong features leave their mark on the 
spectrum. 

Some weak features at about 2.1 and 2.38 /xm are due to Mg II (21369, 21432, 24041, 24044, 
24125 A) and Si II (21920, 21990 A). The strong, broad feature between 1.62 and 1.75 /um is a blend 
produced by Mg II lines (16760, 16799, 174119, 17717 A) and Si II (16907, 16977, 17183, 17184 A). 
Its blue edge is determined by Mg II. A similar feature is also present in normal luminosity SNe la 
but Mg II does not contribute because the corresponding spectra are formed in layers of explosive 
oxygen burning. 

The strong feature that peaks near 1.35 ^m is due to Si II (13650 A), and weaker features 
with peaks near 1.23 jam and 1.43 jam are a Si II multiple! (13692 to 13696 A) and a blend of O I 
(I 4542 A) and Si II (14454 A), respectively. A very weak feature in the synthetic spectrum near 
1.23^m can be attributed to O I (13164 A), but is well below the noise level of the observed spectra. 

Features with P-Cygni absorption minima at 1.12, 1.06, 1.03 and 1.0 ^m are produced by 
Mg II (11620, 11600 A), O I (11302, 11286 A), Mg II (10914, 10915 & 10950 A), C I (10683 k 
10691 A), and Mg II (10092 A), respectively. The multiple! of Mg II 1.09 p.m is also prominent in 
normal luminosity SNe la but the other features are weak and blended with Fe II, Co II and Ni II. 

The observed spectra go down to about 9700 A. However, there are some interesting features 
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at shorter wavelengths which are worth mentioning. Foremost, the model predicts a very strong 
feature at about 9100 A due to C 1 (9405 A). Its velocity provides an important constraint on the 
minimum velocity for the unburned region, and possible mixing. Other features are produced by 
Mg II (9344 A), Ca II (IR-triplet), and O I (8446 A) and Mg II (8246 A). 

Despite the overall agreement, some discrepancies remain. First, the C I line at 1069lA is 
somewhat too strong, and its Doppler shift on May 6 is too large by about 1500 km s“^. This 
feature is formed well above the photosphere (13000 km s“^). This might indicate that SN 1999by 
has a slightly lower temperature in the outer layers compared to the model, which would cause less 
excitation of carbon and, might cut off the high velocity contribution to the absorption. Secondly, 
in the model spectrum there seems to be a lack of emission produced by Si II near 1.18 ^m. We 
note that similar problems with Si are also commonly present in the optical Si II line at about 
6380 A for normal SNe la (Hoflich 1995; Lentz et al. 2000). They may be related to the excitation 
process of Si but their origin is not well understood. 


4-2.2. IR-spectra after maximum light 

The spectrum from May 16 (Fig. 12) is formed in layers of explosive oxygen burning. The 
lines can be attributed to Mg I, Ca II, Si II and S I. All features due to Ca II and Si II seen in the 
pre-maximum spectra can still be identified. No lines due to O I and C I are present, but strong 
lines of Ca II and S I appear. Strong features at about 1.14 //m and 1.1 ^um are due to Ca II lines 
(11784, 11839 & 11849 A) and S I (10821 A) & Mg I (11800 A), and the broad feature at 0.91 ^m 
(not covered in this spectrum, but included for completeness) is a blend of Ca II (9890 A) and 
Mg II (lOOOOA). Weaker lines (Fe II, 9997A; S I, 10455 A) contribute to blends around 1.03 ^m. 

Some weak lines due to Co II and Fe II lines occur throughout the spectrum but the most 
characteristic feature in post maximum IR-spectra of normal SNe la is missing. Starting a few days 
after maximum, normal SNe la show a wide emission feature between 1.5 and 1.8 pLia. produced 
by thousands of Fe II, Co II Sz Ni II lines. Typically, their emission flux is about twice as large 
as the adjoining continuum. This feature marks the lower velocity end of the region of complete 
silicon burning, and it seems to be common in all normal luminosity SNe la such as SN 1986G, 
SN 2000br and SN 2000cx (Wheeler et al. 1998; Marion 2001, priv. comm.; Rudy et al. 2001). In 
5p0z22.8, this broad emission feature is not seen because the photosphere has yet to recede to the 
layers of complete silicon burning. On May 24th (Fig. 13) this feature finally appears but is very 
weak. Overall, two weeks after maximum light, the spectrum in dominated by a large number of 
Fe II and Co II lines. Strong features of Fe II are around 9600 A 1.2 ^m and 1.4 ^um with strong 
Co II lines from 9500 to 12443 A, and 14000A - 16700 A. 

In our synthetic spectrum the Fe II, Co II &: Ni II feature 1.6 to 1.8 lam does not show a local 
minimum at ~ 1.7/rm as observed in SN1999by, and as also predicted by normal bright models for 
SN1986G (Hoflich 1997; Wheeler et al. 1998). In our subluminous model, the emission near 1.7 urn 
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comes from the central Ni-rich region. The absence of such emission in the observed spectrum may 
indicate that the very central region (M(r) < 0.2Mq) of the model is too transparent compared 
SN1999by. This discrepancy is not critical for measuring the chemical structure of the envelope 
because this is determined by the blue edge. The difference in opacity may be caused by missing 
iron lines in the line list (§3.1.3), or an underestimation of the excitation by non-thermal electrons 
because we assume local deposition for the electrons from the /3~^ decay of ^^Co. 


4.3. Mixing processes 

As we have seen, model 5p0z22.8 reproduces the basic features of both optical and IR spectra, 
including the time evolution. This suggests that DD models are a reasonably good description for 
subluminous SNe la. 

Unfortunately, we are still missing a handle on an important piece of physics. The propaga¬ 
tion of a detonation front is well understood but the description of the deflagration front and the 
deflagration to detonation transition (DDT) pose problems. Currently, state of the art allows us 
to follow the front only through the phase of linear instabilities, i.e. only the early part of the de¬ 
flagration phase. The resulting structures of these calculations cannot account for the observations 
of typical SNe la (Khokhlov 2001), because a significant fraction of the WD (~ O.SMq) remains 
unburned. 

Nevertheless, these calculations point toward a more general problem, or toward a better un¬ 
derstanding of the nature of subluminous SNe la. On a microscopic scale, a deflagration propagates 
due to heat conduction by electrons. Though the laminar flame speed in SNe la is well known, 
the front has been found to be Rayleigh-Taylor (RT) unstable, increasing the effective speed of the 
burning front (Nomoto, Sugimoto, & Neo 1976). More recently, significant progress has been made 
toward a better understanding of the physics of flames. Starting from static WDs, hydrodynamic 
calculations of the deflagration fronts have been performed in 2-D (Niemeyer &: Hillebrandt 1995; 
Lisewski et al. 2000), and full 3-D (Khokhlov 1995, 2001). 

The calculations by Khokhlov (2001) demonstrated a complicated morphology for the burning 
front. Plumes of burned material will fill a significant fraction of the WD, and unburned or partially- 
burned material can be seen near the center. Thus iron-rich elements will not be confined to the 
central region as in 1-D models. Khokhlov (2001) finds that while the expansion of the envelope 
becomes almost spherical, the inhomogeneous chemical structure will fill about 50 to 70% of the 
star (in mass). If a DDT occurs at densities needed to reproduce normal-bright SNe la, most of 
the unburned fuel in these regions will be burned to iron-group elements during the detonation 
phase. This will eliminate the chemical inhomogeneities. However, in our subluminous model, the 
chemical structure imposed during the deflagration phase must be expected to survive because no 
or very little is produced during the detonation. 

In order to test for this effect, we have mixed the inner layers of model 5p0z22.8 up to an 
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expansion velocity of 8000 km s“^. This brings iron-group elements into the outer regions and 
mixes intermediate mass elements (Si, S, Ca and some O) into the inner, hotter layers. The layers 
with expansion velocities less than 8000 km s“^ become visible a few days after maximum light 
(see above). In Figure 14, a comparison is shown between the spectrum of SN1999by on May 16th, 
1999, and the synthetic spectrum about one week after maximum light, and problems with the 
spectral fit are apparent. As can be expected from the discussion above, the effective extension of 
the iron-rich layers to 8000 km s“^ causes the 1.5 to 1.8 /xm blend to become rather strong, and 
line blanketing due to iron-group elements (Fe II Sz Co II) become too strong in the region around 1 
to 1.2 ^m. Moreover, at about 0.99, 1.02, 1.05 & 1.16 /um, (absorption components) strong blends 
can be seen due to O I (11302, 11287 & 11297A), Si II (11715, 11745, 13650 k 13650 A), Ca II 
(9890, 9931, 9997, 11839 k 11745 A) and S I (10455 k 108212 A). 

We have to take into account that our approach assumes microscopic mixing whereas RT 
instabilities provide large scale inhomogeneities. If and the intermediate mass elements are 
separated, direct excitation of intermediate mass elements may be reduced. With respect to the 
radiation transport, the main difference is the covering factor. This will effect line blanketing and 
thus the line blanketing between 1 to 1.2 ^m may not pose a problem. However, the emission 
features between 1.5 to 1.8 ;um will not go away. The reduced excitation by non-thermal electrons 
and y’s in case of macroscopic mixing are important at low optical depths. Less excitation may 
reduce the problem with the O, Si, S & Ca lines but the problem will no disappear because the 
corresponding features are formed close to the photosphere where the excitation is thermal in 
nature. 

These problems strongly suggest that large scale mixing did not occur in SN 1999by. This may 
provide a key for our understanding of the nature of these objects. 


5. Implications for the SN la progenitors, and alternative scenarios 

As previously discussed, the pre-conditioning of the WD before the explosion may provide a 
key for understanding the nature of subluminous SNe la. Such pre-conditioning may include the 
main sequence mass of the progenitor mass and its metallicity (Dominguez, Hoflich, k Straniero 
2001; Hoflich, Wheeler, k Thielemann 1998; Iwamoto et al. 1999), the accretion history (Danger 
et al. 2000), or large-scale velocity fields such as turbulence prior to the runaway (Hoflich k Stein 
2001). Up to now, none of these effects have been included in detailed 3-D calculations. Therefore, 
the suggestions (below) require further investigation before taken as explanation. 

Beyond its subluminous nature, SN 1999by showed another peculiarity with respect to normal- 
bright SNe la. The observed polarization of SN 1999by was rather high (P ~ 0.7%) with a well- 
defined axis of symmetry, whereas normal SNe la show little or no polarization (P < 0.2%; Wang, 
Howell, Wheeler, k Hoflich 1997). ^From a detailed analysis of the polarization spectra, Howell 
et al. (2001) concluded that the overall geometry of SN 1999by showed a large scale, probably 



- 19 - 


rotational asymmetry of about 20%. The authors suggested several explanations. Either, SN 1999by 
may be the result of the explosion of a rapidly rotating WD close to the break-up, or a result of 
a merger scenario. In a rapidly rotating WD, large-scale circulations (Erigichi & Muller 1985) 
might influence the deflagration phase of explosive burning, perhaps by breaking up large eddies 
and preventing their rise. 

Alternatively, it should be remembered that the chemical structure of the exploding WD 
depends on the pre-expansion of the WD prior to the detonation phase. More precisely, we need 
to burn about 0.2Mq to lift the WD in its gravitational potential. This mass is significantly more 
than is burned during the smoldering phase (the phase of slow convective burning just prior to the 
explosion). Even in the presence of strong turbulence, the amount of burned material is much less 
(< 1O“^M0, Hoflich & Stein 2001). In all models for the ignition (e.g. Nomoto 1982; Garcia- 
Senz &: Woosley 1995), the central temperature rises by non-explosive carbon burning during the 
last few minutes to hours before the runaway because the heat cannot be dispersed. It may be 
feasible that large scale velocity fields due to rotation may extend this phase, and significantly 
increase the amount of burning prior to the explosion. Thus, the pre-expansion in subluminous 
SNe la may occur directly during the smoldering phase, and the deflagration phase may be even 
skipped, i.e. a detonation occurs promptly, and we may have a smoldering detonation model (SD). 
In SD-models, no is produced during the phase of pre-expansion and, consequently, this 

mechanism does not predict a lower limit for and, thus. My. In principle, SD-models may be 
an attractive alternative also for normal-bright SNe because they omit the need for a DDT but, 
from the present understanding of the runaway, in general, extensive burning prior to the runaway 
cannot be expected (see above). From the point of the spectral analysis, this option is intriguing 
because burning during the smoldering phase precedes only to O and the time scales of burning 
are sufficient to completely homogenize the chemical structure of the WD. 

Finally, there is the question of alternative scenarios for SNe la in the context of subluminous 
SNe la. Pure deflagrations of Mc</j-WDs without any DDT may produce little but, according 
to current simulations, a large fraction of the WD would remain unburned. Even if this problem 
were solved in the future, the problem of inhomogeneous chemistry and Si-rich layers close to 
the center would remain. Helium triggered detonations of sub-Chandrasekhar mass WDs have 
been suggested in the past. However, their distinct feature is an outer layer of above a low 
mass layer of explosive oxygen burning, a layer of incomplete Si burning, and some close to 
the center. This chemical structure is in strong contradiction to the optical spectra and LCs of 
SN 1999by, and their evolution with time. 

Alternatively, setting aside the problem of triggering the explosion, central detonations in low- 
mass WDs have been suggested in the past (e.g., Ruiz-Lapuente et al. 1993). To produce little 
^^Ni, the explosion must occur in WDs with central densities < 1 x 10^ g cm“^ which corresponds 
to M{WD) < 0.7Mq. In contrast to the DD-models, these models would not show a lower limit for 
the production and My, and they would avoid the problem of mixing during a deflagration 

phase. Nonetheless, we regard these models as unlikely candidates for SN 1999by because there are 
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a couple of severe problems beside the lack of a triggering mechanism for the explosion. The total 
mass of the ejecta is lower (by a factor of 2) than is required to keep the line forming region in the 
layers of explosive carbon and incomplete Si-burning up to two and four weeks after the explosion, 
respectively. The specific production of nuclear energy does not depend on Mwr,. In this case, the 
diffusion time scales and, thus, the rise times to maximum light are oc M\yd (e-g- Pinto & Eastman 
2000). Compared to the LCs for a WD with Mch, we expect rise times to be shorter by ~ 2 which 
is incompatible with SN 1999by observations. 

This leaves the merging scenario as an alternative as we may expect a significant asphericity in 
the explosion if the merging occurs on hydrodynamical time scales. Nevertheless, this alternative is 
not a very attractive one. Based on 1-D parameterized models, Khokhlov, Muller, & Hoflich (1993) 
provided chemical structures for explosions which may resemble merger models. As a general 
feature, these models show a shell-like structure with a massive outer layer of unburned C/0 which 
encloses a thin, low mass layer of explosive carbon burning. This pattern is not consistent with 
the IR-spectra for SNe 1999by up to maximum light. However, it should be kept in mind that, 
besides the limitations of 1-D models for mergers, no chemical structures are available for merging 
subluminous models. A more severe problem for mergers may be that merger models tend to ignite 
prematurely and produce a detonation wave that burns the C/0 WD to a 0/Ne/Mg WD prior to 
the completion of the merger. Thus, merging may result in an accretion induced collapse rather 
than a SNe la (Saio & Nomoto 1998). 


6. Summary 

IR-spectra of the subluminous SN1999by have been presented which cover the time evolution 
from about 4 days before to 2 weeks before maximum light. This is the first subluminous SN la 
(and arguably the first SN la) for which IR spectra with good time coverage are available. These 
observations allowed us to determine the chemical structure of the SN envelope. 

Based on a delayed detonation model, a self-consistent set of hydrodynamic explosions, light 
curves, and synthetic spectra have been calculated. This analysis has only two free parameters: 
the initial structure of the progenitor and the description of the nuclear burning front. The light 
curves and spectral evolution follow directly from the explosion model without any further tuning, 
thus providing a tight link between the model physics and the predicted observables. 

By varying a single parameter, the transition density at which detonation occurs, a set of 
models has been constructed which spans the observed brightness variation of Type la supernovae. 
The absolute maximum brightness depends primarily on the production which, for DD-models, 
depends mainly on the transition density ptr- The brightness-decline relation My{/S.M^t=ibd) 
observed in SNe la is also reproduced in these models. In the DD scenario this relation is a result 
of the temperature dependence of the opacity, or more precisely, as a consequence of the rapid drop 
in the opacities for temperatures less than about 7000 to 8000 K. 
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Within Mch WD models, a certain amount of burning during the deflagration phase is needed 
to pre-expand the WD and avoid burning the entire star to ^^Ni. This implies a lower limit for 
the mass of about O.IMq ±25% and, consequently, implies a minimum brightness for SNe la 
within this scenario. The best model for SN 1999by {ptr = 8 x 10^ g cm“^, selected by matching 
to the predicted and observed optical light curves) is close to this minimum Ni yield. The data are 
consistent with little or no interstellar reddening {E[B — V) < 0.12™), and the derived distance 
is 11 ± 2.5 Mpc or 12 ± 1 Mpc if we take the limit for E(B — V) from the models or assume 
E{B — V) = 0.015™ according Schlegel et al. (1998) for the galaxy. 

Without any further modification, this subluminous model has been used to analyse the IR- 
spectra from May 6, May 10, May 16 and May 24, 1999, which correspond to —4 d, 0 d, ±7 d and 
±14 d after maximum light. The photosphere {rThomson = 1) recedes from 14000, to 10500, 6500 
and 4000 km s“^. The observed and theoretical spectra agree reasonably well with respect to the 
Doppler shift of lines and all strong features could be identified. 

Before maximum light, the spectra are dominated by products of about explosive carbon 
burning (O, Mg), and Si. Spectra taken at ±7 d and ±14 d after maximum are dominated by 
products of incomplete Si burning. At about 2 weeks after maximum, the iron-group elements 
begin to show up. The long duration of the phases dominated by layers of explosive carbon burning 
and incomplete Si burning implies massive layers of these burning stages that are comparable with 
our model results. (~ 0.45 and 0.65Mq for explosive carbon and incomplete silicon burning, 
respectively.) This, together with the mass, argues that SN 1999by was the explosion of a 
WD at or near the Chandrasekhar mass. 

Finally, we note that the observed IR spectra are at odds with recent 3-D calculations (Khokhlov 
2001) which predict that large scale chemical inhomogeneities filling 50 - 70% (in mass) of the WD 
will be formed during the deflagration phase. When the effect of such inhomogeneous mixing is 
tested by mixing the inner layers of our SN 1999by model, significant differences appear between 
the model spectra and the observed data. This suggests that no significant large-scale mixing took 
place in SN 1999by. The lack of observed mixing and the asphericity seen in SN 1999by may be 
important clues into the nature of subluminous SNe la, and may be related to the reason for a low 
DDT transition density (and hence the low luminosity). Alternatively, we may have an extended 
smoldering phase of the WD prior to the explosion which skips the deflagration phase altogether 
(see §5). In either case, the pre-conditioning of the WD prior to the explosion seems to be a key 
for understanding SNe la. 
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Table 1. Log of Near-Infrared Observations 


Date 

Epoch 

Setup 

Exp. Time 

Observer 

6 May 

-4 d 

0.96-1.8 ;um 

4x600 s 

Gerardy 



1.95-2.4 

3x600 s 

5? 

7 May 

-3 d 

0.96-1.8 fiui 

3x600 s 

5? 



1.95-2.4 fiui 

3x600 s 

5? 

8 May 

-2 d 

0.96-1.8 fiui 

3x600 s 

5? 

9 May 

-1 d 

0.96-1.8 

3x600 s 

5? 



1.2-2.2 fim. 

3x600 s 

5? 

10 May 

0 d 

0.96-1.8 jj-m 

3x600 s 

5? 

D 

5? 

1.2-2.2 |Um 

3x600 s 

5? 

5? 

?? 

1.95-2.4 fiui 

3x600 s 

5? 

16 May 

6 d 

0.96-1.8 jj-m 

3x600 s 

Martini & Steed 

18 May 

8 d 

0.96-1.8 jj-m 

3x600 s 

5? 

20 May 

10 d 

0.96-1.8 /um 

3x600 s 

5? 

21 May 

11 d 

0.96-1.8 jj-m 

3x600 s 

Sakai 

24 May 

14 d 

0.96-1.8 jj-m 

3x600 s 

5? 

NIR Photometry 

Date 

Epoch 

J 

H 

K 

6 May 

-4 d 

13.28 ±0.09 

13.23 ±0.05 

13.19 ±0.05 

7 May 

-3 d 

13.28 ±0.09 

13.21 ±0.05 

13.14 ±0.05 

8 May 

-2 d 

13.07 ±0.09 

13.04 ± 0.05 











Table 2. 


Properties of the models 


Name 

ptr[lQ^g/cm 3] 

[Mq] 

Mv 

tv 

B-V 

AMi5 

5p0z22.8 

8. 

0.095 

-17.21 

13.9 

0.47 

1.45 

5p0z22.10 

10. 

0.107 

-17.35 

14.1 

0.42 

1.37 

5p0z22.12 

12. 

0.137 

-17.63 

14.7 

0.38 

1.32 

5p0z22.14 

14. 

0.153 

-17.72 

14.9 

0.22 

1.30 

5p0z22.16 

16. 

0.268 

-18.72 

15.8 

0.14 

1.26 

5p0z22.18 

18. 

0.365 

-18.82 

16.6 

0.08 

1.21 

5p0z22.20 

20. 

0.454 

-18.96 

17.0 

0.02 

1.19 

5p0z22.23 

23. 

0.561 

-19.21 

18.2 

-0.02 

1.05 

5p0z22.25 

25. 

0.602 

-19.29 

18.6 

-0.02 

1.00 

5p0z22.27 

27. 

0.629 

-19.35 

19.0 

-0.03 

0.91 
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Fig. 1.— Near-infrared spectral evolution of SN1999by. Epochs are relative to the date of Fiviax, 
May 10, 1999. For clairity, the spectra have been shifted vertically, and regions of very low S/N 
due to telluric absorption have been omitted. 




p /p(central) p /p(central) p /p(central) 


~ 30 - 





o 

<u 

CO 


B 


CO 

O 


> 



O 

<D 

CO 


B 


CO 

O 


> 


Fig. 2.— Density (blue, dotted) and velocity (red, solid) as a function of the mass for models with 
ptr = 8, 16 and 25 x 10®(from top to bottom). 
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5p02822 with = 16c E6g/ccm [Stable] 



vel [1000 km/sec] 



Fig. 3.— Abundances of stable isotopes as a function of the expansion velocity for models with 
ptr = 8, 16 and 25 xlO®g/cm“^ (from top to bottom). In addition, is given. The curves with 
the highest abundance close to the center correspond to ^^Fe, and ^^Fe. 
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Fig. 4.— B (left) and V (right) light curves for models with ptr = 8, 12, 16, 20 and 25 x 10^g/cm^ 
(from bottom to top). 
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Fig. 5.— Maximum brightness My as a function of the mass (upper left), ptr (lower left), 

and Mi/(AMAt=i 5 d) (lower right), and {B — V) as a function of My (upper right). For all panels 
except the lower right, the points correspond to models with ptr of 8, 10, 12, 14, 16, 18, 20, 23, 25 
and 27 xlO^g/cm^ from left to right. In the lower right panel, this order is reversed. 
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Fig. 6 .— Temperature (left), Rosseland optical depth (center), and energy deposition by 7 -rays 
(right), as functions of radius at day 10, 14, 21 and 28 (from top to bottom). 
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Spectral Evolution of the Model 



at 10.6, 15, 22 and 29 days after the explosion. The flux has been normalized at 7000 A. 
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Fig. 8.— Comparison of the observed B (left) and V (right) light curves of SN1999by with the 
predicted light curves of model 5p0z22.08. 
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SN1999by on 10may99 



wavelength[A] 

Fig. 9.— Comparison of the optical spectrum of SN1999by on May 10, 1999, (blue) with the 
theoretical spectrum of 5p0z22.8 at maximum light (red, 15 days after the explosion). 
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SNig99by on 06may99 




theoretical IR-spectrum of 5p0z22.8 ~ 4 days before maximum light (10 days after the explosion). 
Between 1.35 and 1.45 ;um, and longward of 2.45 /rm, the S/N in the observed spectrum is very 
low due to strong telluric absorption. 




























- 39 - 


SN1999by on 10may99 



SN1999by on 10inay99 



theoretical IR-spectrum of 5p0z22.8 at maximum light (15 days after the explosion). Between 1.35 
and 1.45 pm, between 1.8 and 1.9 pm, and longward of 2.45 pm, the S/N in the observed spectrum 
is very low due to strong telluric absorption. 
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SNI999by on 16may99 



Fig. 12.— Comparison of the observed NIR spectrum of SN1999by on May 16, 1999, with the 
theoretical IR-spectrum of 5p0z22.8 seven days after maximum light (22 days after the explosion). 
Between 1.35 and 1.45 /rm the S/N in the observed spectrum is very low due to strong telluric 
absorption. 
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SN1999by on 24may99 



wave]ength[A] 

Fig. 13.— Comparison of the observed NIR spectrum of SN1999by on May 24, 1999, with the 
theoretical IR-spectrum of 5p0z22.8 14 days after maximum light (29 days after the explosion). 
Between 1.35 and 1.45 /rm the S/N in the observed spectrum is very low due to strong telluric 
absorption. 
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SN1999by on 16inay99 



wavelengtli[A] 


Fig. 14.— Comparison of the observed NIR spectrum of SN1999by on May 16, 1999, with the 
theoretical IR-spectrum of 5p0z22.8 seven days after maximum light if we impose mixing in the 
inner 0.7 Mq. Line blanketing due to a large number of Fell and Coll lines becomes strong as on 
May 24, 1999 in the unmixed model (see Fig. 12). 
















































